LPT Orsay 99/19 



hcp-ph/990333S 



(N 



Event— by— event fluctuations in heavy— ion collisions and the 

quark— gluon string model 

A. Capella, 1 , E. G. Ferreiro 1 , A. B. Kaidalov 2 , 

1 Laboratoire de Physique Theorique, 
Unite Mixte de Recherche (CNRS) UMR 8627, 
Universite de Paris XI, Bat. 210, 91405 Orsay Cedex, France, 



> 
OO . 

£f) . Institute of Theoretical and Experimental Physics, 

m _ 

^ . B. Cheremushinskaya 25, 111 259 Moscow, Russia 

ON 



Abstract 



X 

d ' We apply dynamical string models of heavy-ions collisions at high energies 

to the analysis of event-by-event fluctuations. Main attention is devoted to a 
new variable proposed to study "equilibration" in heavy-ions collisions. Recent 
results of the NA49 collaboration at CERN SPS are compared with predictions 
of the Quark-Gluon String Model (QGSM), which gives a good description of 
different aspects of multiparticle production for collisions of nucleons and nuclei. 
It is shown that the new observable and other results of the NA49 analysis of 
event-by-event fluctuations are correctly reproduced in the model. We discuss 
dynamical effects responsible for these fluctuations and give the predictions for 
p-p, p-Pb and Pb-Pb collisions at RHIC and higher energies. 
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I. INTRODUCTION 



An event-by-event analysis of heavy-ions collisions can give important information on the 
dynamics of these processes. In the paper |] it was proposed to use an event-by-event analysis 
of transverse momentum fluctuations as a method for the study of "equilibration" in high- 
energy nucleus-nucleus collisions. For this purpose a special variable $ has been introduced 
in ref. |l| (the definition of this variable will be given below). The problem of "equilibration" 
in high-energy heavy-ions collisions is very important in order to understand the conditions 
for quark-gluon plasma formation. Recent experimental results on event-by-event analysis 
of Pb-Pb collisions at CERN SPS by the NA49 collaboration show that the value of $ is 
substantially smaller than expected in the case of independent nucleon-nucleon collisions and 
its smallness was considered as an indication of "equilibration" in the system. 

This result has been discussed in the framework of different theoretical models. It was 
shown || that the increase of transverse momenta of hadrons due to multiple rescatterings 
leads to a substantial increase of $. Incorporating this effect in the model of ref. Q one finds 
an even stronger disagreement with the NA49 result. It was also demonstrated that strings 
fusion also leads to an increase of $ in disagreement with experiment The results on the 

influence of final states interactions on the observable $ are contradictory: in ref. || it was 
shown that final states interactions in the framework of the string model of ref. [Bj have a small 
effect and do not allow to reach agreement with experiment, while in ref. |7j it was argued that 
final states interactions in the framework of the UrQMD model are essential and can decrease 
$ to a value consistent with experiment. On the other hand, it was shown in ref. || that, in the 
case of fully equilibrated hadronic gas made mostly of pions, one expects large positive values 
for the variable $ not consistent with the experimental observation. 

In this note we study event-by-event fluctuations using the Monte-Carlo formulation |§ 
of the Quark-Gluon String Model (QGSM) . The QGSM and the Dual Parton Model 
(DPM) |ll|) are closely related dynamical models based on 1/N-expansion in QCD, string 
fragmentation and reggeon calculus. They give a good description of many characteristics of 

2 



multiparticle production in hadron-hadron, hadron-nucleus and nucleus-nucleus collisions (for 
a review see refs. ||12||). Nuclear interactions in this model are treated in the Glauber-Gribov 
approach. It will be shown that the model reproduces the results of the event-by-event analysis 
of the NA49 experiment for the quantity $ as well as for other fluctuations observed in this 
experiment 0. We analyze the reason for the decrease of the quantity $ from p-p to Pb-Pb 
collisions seen by the NA49 experiment and come to the conclusion that the quantity $ is 
sensitive to many details of the interaction and can hardly be considered as a good measure 
of "equilibration" in the system. We predict a strong increase of the value of $ at energies of 
RHIC and higher. The model also gives definite predictions for event-by-event fluctuations in 
p-Pb collisions. 

II. ANALYSIS OF EVENT BY EVENT FLUCTUATIONS OF TRANSVERSE 

MOMENTA 

Let us remind briefly the method to study event-by-event fluctuations of transverse mo- 
menta of produced particles introduced in ref. It was proposed to define for each particle in 
a given event a variable Zi = pn— < Pt >, where pxi is the transverse momentum of the parti- 
cle i and < > is the mean transverse momentum of particles averaged over all events. Using 
Zi the quantity Z = Yif =l Zi is defined, where N is the total number of particles in the event. If 
nucleus-nucleus collisions can be considered as a superposition of independent nucleon-nucleon 
collisions then it can be shown [jl| that 



< N > AA < N > NN 

A derivation of this result is given in Appendix 1. 

The averaging in eq.(l) is over all events in a given kinematical region. It was proposed 
in ref. to characterize the degree of fluctuations by the variable 



< Z 2 > AA < Z 2 



>NN 



(1) 




(2) 
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where < z 2 > is the second moment of the single particle inclusive z-distribution. The quantity 
< z 2 > corresponds to purely statistical fluctuations and is determined by mixing particles 
from different events. It was emphasized in ref. that, if nucleus-nucleus collisions would be 
a simple superposition of independent nucleon-nucleon collisions, then the variable $ would be 
the same as in the nucleon-nucleon case (see Appendix 1). In nucleon-nucleon collisions the 
quantity $ is different from zero due to dynamical correlations and, in particular, due to the 
dependence of < > on the number of produced particles. It was proposed to attribute a 
possible decrease of the quantity $ in A-A collisions to the effects of "equilibration". 

Let us note that the model of independent N-N collisions for nucleus-nucleus interactions 
is an extremely oversimplified one. The Glauber model at high energies is not equivalent to in- 
dependent N-N collisions even for N-A interactions. The space-time picture of hadron-nucleus 
interactions at high energies is absolutely different from a simple picture of successive reinter- 
actions of an initial hadron with nucleons of the nucleus (see e.g. |13|-p^|). For nucleus-nucleus 



interactions there are extra correlations [14,16]. The model of independent N-N collisions does 



not even satisfy energy-momentum conservation, as a nucleon of one nucleus can not interact 
inelastically several times with nucleons of another nucleus having the same energy at each 
interaction. 

In the QGSM as well as in DPM, the effects of multiple interactions in hadron-nucleus 
and nucleus -nucleus collisions are taken into account in the approach based on the topological 
expansion in QCD [|l^]. Probabilities of rescatterings are calculated in the framework of the 
Glauber-Gribov theory and multiparticle configurations in the final state are determined using 
AGK [|17] cutting rules. In these models the Pomeron is related to the cylinder type diagrams, 
which correspond to the production of two chains of particles due to decays of two qq — q strings. 
Multi-Pomeron exchanges are related to multi-cylinder diagrams which which produce extra 
chains of type q — q- They are especially important in interactions with nuclei. Fragmentation 



of strings into hadrons is described according to "regge counting rules" |18j, which give correct 
triple-regge and double-regge limits of inclusive cross sections. Note, however, that in the 
Monte Carlo version used in this paper, the Artru-Menessier string fragmentation scheme is 



implemented instead. All conservation laws (including energy-momentum conservations) are 
satisfied in this approach. 

Let us emphasize that at energies y/s ~ 10 GeV the cylinder-type diagrams give the 
dominant contributions for N-N collisions. Extra q — q chains due to multi-cylinder diagrams 
have rather small length in rapidity (short chains) and do not lead to substantial contributions 
to particle production. In nucleon-nucleus and nucleus-nucleus collisions, the number of short 
chains is strongly increased compared to the nucleon-nucleon case (it is proportional to a 
number of collisions) and they should be taken into account in any realistic calculations of 
multiparticle production on nuclei. This means that for p-A and A-B collisions there are extra 
"clusters" of particles (short chains, of type q — q) compared to the nucleon-nucleon interactions 
"clusters" (long chains connecting valence quarks and diquarks of the colliding nucleons). 

So we come to the conclusion that in the relativistic Glauber-Gribov dynamics the char- 
acteristics of final particles in N-A and A-B collisions can not be expressed in terms of N-N 
collisions only, as it was assumed in ref. and eq.(l) is not valid in general. In Appendix 1 we 
give as an illustrative example the results in a model with two types of clusters. This model is 
a generalization of the single cluster model of ref. 0, and is much closer to QGSM and DPM. 

The results of the Monte Carlo calculation for the quantity $ are shown in Table 1 for p-p, 
p-Pb and central Pb-Pb collisions at SPS energies (y/s = 19.4 GeV) and at RHIC (y/s = 200 
GeV). Predictions of the model for $ are quite different for these two energies. At SPS there 
is a strong reduction of the quantity $ for nuclear collisions compared to N-N collisions, while 
at RHIC the quantity <3> is predicted to be much larger than at SPS and about the same for 
Pb-Pb and p-p collisions. At LHC energies the value of obtained in our model, is 160 MeV 
for p-p and even larger for Pb-Pb collisions (the Monte Carlo code we use does not allow to 
calculate precise values of the correlations at LHC energies due to a too large number of particle 
produced in each event; so in Table 1 we give predictions of the model at y/s = 540 GeV and 
y/s = 1 TeV to show the energy dependence of fluctuations). The results for SPS energies are 
in a reasonable agreement with experimental data of the NA49 Collaboration ||. The result 
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for p-p interactions at this energy is even higher than the estimate, based on the dependence 
of < pt > on the number of charged particles, given in ref. |]. The model reproduces this 
correlation reasonably well and shows that this is not the only source of fluctuations leading to a 
non-zero value of $. We find that the quantity <3> is sensitive also to other types of correlations 
and, in particular, to the correlations related to conservation of px in the process. 

Let us note that the quantity $ at SPS energies is very small and is defined in eq.(2) 
as a difference of two large numbers (see Table 1), so it is very sensitive to all details of 
dynamical models. Because of its smallness it is difficult to obtain a good accuracy in a Monte 
Carlo calculation of this quantity (especially for nucleus-nucleus collisions, where the maximum 
statistics possible in the Monte-Carlo is ~ 5000 events). In order to increase the statistics and 
to reduce this uncertainty we give in Table 1 results obtained for the total rapidity interval, 
while the experimental data of the NA49 Collaboration were obtained in a fixed rapidity interval 
4 < y l ® < 5.5. The error in the values of $ in Table 1 is about 1 MeV for the lowest energies, 
and increases at high energy. Other properties of event by event fluctuations observed by NA49 
Collaboration [[| were calculated under the condition of the experiment and are reproduced by 
the model reasonably well as it is shown in Fig.l. 



from p-p to Pb-Pb collisions, but there is an even larger increase for V < z 2 > due to the 
increase of < px > and to a change in the form of the pt distribution. The effect of the 
correlations between < px > and the number of charged particles due to rescatterings is, to a 
large extent, compensated at these energies by energy-momentum conservation effects. As a 
result, we find no dependence of < px > on n c h for Pb-Pb collisions at SPS. For RHIC energies 
a strong increase in the values of both yj and especially of <3> is predicted (see Table 1). At 
these energies the increase of average transverse momentum with the number of rescatterings 
becomes very important in p-p interactions and is reproduced by the QGSM (Fig. 2a). It is 
shown in Fig. 2b that an increase of < px > with multiplicity is predicted at these energies 
even for Pb-Pb collisions, although the effect is less pronounced for heavy-ions collisions than 
for p-p. At LHC energies all these effects will be stronger than at RHIC and will produce an 



It follows from Table 1 that at SPS energies there is an increase of the quantity 
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increase of the quantity $ in p-p, p-A and A-A collisions as energy increases. 

The predictions of the model can be easily tested in future experiments at RHIC and 

LHC. 

III. CONCLUSIONS 

We have shown, in a Monte Carlo version of the QGSM, that at SPS energies the quantity 
$, caracterizing event-by-event transverse momentum fluctuations, decreases from $ ~ 9 MeV 
in p-p collisions to $ ~ 2 MeV in central Pb-Pb collisions. This result for Pb-Pb collisions 
agrees with the measurement of the NA49 Collaboration. In ref. 0, such a decrease between 
p-p and Pb-Pb was considered to be a test of equilibration of the dense system produced 
in central heavy ion collisions. We have obtained the same result in the framework of an 
independent string model. 

At RHIC energies, we predict an increase in the value of $ ($ = 75 4- 80 MeV). In this 
case $ will be approximately the same in p-p and central Pb-Pb collisions. At higher energies 
the value of $ is predicted to be larger and to increase from p-p to central Pb-Pb collisions. 

Our analysis indicates that the quantity $ can hardly be considered as a good measure of 
"equilibration" in the system. However, it can be used as a sensitive test of dynamical models. 
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APPENDIX 1 



Here we will consider a simplified model of multiparticle production with two types of 
"clusters", which is a generalization of the model of ref. where clusters of only a single 
type were produced. As discussed above, these " clusters" correspond to qq — q chains (clusters 
of the first type) and q — q chains (clusters of the second type). Nucleon-nucleon collision 
at SPS-energies can be described with a good accuracy by the production of two clusters of 
the first type, while for proton-nucleus and nucleus-nucleus interactions, production of the 
second type of clusters is important even in this energy range. For independent production 
of k± clusters of the first type and k 2 - of the second type, with average transverse momenta 
< pt >i and multiplicity < n >i for the cluster i, the following results for < Z > and < Z 2 > 
can be obtained: 

< Z > klk2 = h < Z >i +k 2 < Z > 2 (A.l) 
where < Z >;=< n >< (< p T > t - < P T >) and < P T >= «fci><"i><PT>i+<fc 2 ><n> 2 <p T > 2 ) 

1 1 ' 1 (<fci><n>i + <K2><n>2) 

< Z 2 > klk2 = h < Z 2 >! +k 2 < Z 2 > 2 +k 1 (k 1 - 1) < Z >j + 

+ k 2 (k 2 - 1) < Z > 2 2 +2k t k 2 < Z > x < Z > 2 (A.2) 

The expression for < Z 2 > kl k 2 can be rewritten as 

< Z 2 > klk2 = h(< Z 2 > X -<Z >{) + k 2 (< Z 2 > 2 -<Z > 2 2 )+ < Z > 2 kik2 . (A.3) 

For the quantity $ in eq.(||D it is important that contrary to the case of a single cluster, 
the expression for < Z 2 > contains negative terms proportional to < Z >?. 

Next we take the average over the number of produced clusters with some distribution 

p 
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« Z 2 >k 1 k 2 >— X! Pkik 2 < Z 2 >feifc 2 — 

kl,k2 

--< h > (< Z 2 >! - < Z >j)+ < k 2 > (< Z 2 > 2 - < Z > 2 2 )+ « Z > 2 klk2 > . (A.4) 



In the following we will denote this averaging simply by < Z 2 >. We shall concentrate on 
p-A collisions. In this case, it is easy to show that the last term in eq. ( [A.4] ) is small and can 



be neglected. To prove this we note that for p-A collisions, k 2 = k% — 2 |T2|] with k\ — v + 1, 
where V is the average number of collisions and, thus, 

« Z > 2 klk2 > -«Z > klk2 > 2 = (< k\ > - < h > 2 )(< Z >! + < Z > 2 ) 2 . (A.5) 

Taking into account that (for fixed impact parameter) the distribution in k\ is of a Poisson 
type with (< k\ > — < k\ > 2 ) = c\ < k\ > and that 

<< Z > klk2 >=< ki >< Z >i + < k 2 >< Z > 2 =< h >< Z >! +(< ki > -2) < Z > 2 = we 
obtain: 

« Z >U>= l£L4£Zl . (A . 6) 

felfe 2 < ^ > \ 1 

For large values of < k\ > this quantity is much smaller than the other terms in the right-hand 
side of eq.([PD. 



The expressions for the quantities < N > and < z 2 >, that enter into the definition of $ 
are selfevident 

< >=< ki >< ni > + < k 2 >< n 2 > (A. 7) 



and 

2 < k x >< n >x< z 2 >x + < k 2 X n > 2 < z 2 > 2 

< z A >= . (A. 8) 

< h >< n >i + < k 2 >< n> 2 
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Let us denote 0^ as < z 2 >; (1 + and ^f- = 7*, <^ 2 >< n > 2 = a. Taking into 

<n>i 1 \ v <Z^>; <«i><n>i & 

account that <5, ~ 2-%= and 73 are much smaller than unity we obtain the following approximate 
expression for <3>: 

= [gl ~ 71 ) < ^ >1 + (h - 72)« < z 2 > 2 ] 

where A = (1 + a){< z 2 > x +a < z 2 > 2 ). 

It is important that the terms proportional to 7* give negative contributions to $ and can 
substantially decrease the value of $. 

For the case of clusters of the same type (7, = 0, < Z 2 >i=< Z 2 >2,S\ = 62) we obtain: 

$ = ^V< 2 2 > (A.10) 

both for p-p and p-A. We recover in this way the result of ref. ]!]■ The discussion in this 
Appendix has been restricted to p-A interactions. The situation is more complicated in A-B 
collisions. Actually, even for p-A, we do not claim that the effect discussed in this Appendix 
is the main reason for the decrease of $, obtained in the Monte Carlo calculations (see Table 
1), between p-p and p-A collisions at SPS energies. Nevertheless, our example illustrates the 
important effect that a modification of the model (i. e., going from one to two types of clusters) 
can have on the quantity 
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TABLE CAPTIONS 



Table 1. The results of the Monte-Carlo calculation for the quantities ^|g(MeV), V<^> 
(MeV) and $ (MeV) for p-p, p-Pb and central Pb-Pb-collisions at SPS (y/s = 19.4 GeV), RHIC 
(y/s = 200 GeV) and higher (y/s = 540 GeV and 1 TeV) energies. The results at 1 TeV have 
only an indicative value (see main text). 
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FIGURE CAPTIONS 



Figure 1. Event spectra characterising the multiplicity, transverse momentum and rapidity 
distribution of charged particles per event for Pb-Pb collisions at Pi a b = 158 AGeV/c and 
b < 3.5 fm in the rapidity interval 4 < y < 5.5. The full lines are the Monte Carlo results. 
Experimental data are from ref. |@. 

Figure 2. The dependence of the average transverse momentum on the multiplicity of charged 
particles in the window \rj\ < 2.5 at ^fs = 200 GeV for p— p collisions compared to experimental 
data [|l!J (Fig. 2a) and for Pb-Pb central collisions (Fig. 2b). 
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Table 1 



y/s (GeV) y^^PF ( MeV ) V< z 2 > (MeV) $ (MeV) 



p-p 


19.4 


244.5 


235.5 


9.0 


p-Pb 


19.4 


243.5 


243.0 


0.5 


Pb-Pb 


19.4 


265.6 


263.2 


2.4 


P-P 


200 


387.0 


310.6 


76.4 


p-Pb 


200 


433.7 


367.8 


65.9 


Pb-Pb 


200 


508.9 


429.4 


79.5 


p-p 


540 


450.1 


323.6 


126.5 


p-Pb 


540 


524.3 


397.7 


126.6 


Pb-Pb 


540 


622.6 


475.2 


147.4 


p-p 


1000 


455.5 


324.4 


131 


p-Pb 


1000 


524.5 


397.5 


127 


Pb-Pb 


1000 


704.2 


484.3 


220 
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Figure 1 
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Figure 2 
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